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We calculate the electrical conductivity in the early universe at temperatures below as 
well as above the electroweak vacuum scale, T c ~ lOOGeV. Debye and dynamical screen- 
ing of electric and magnetic interactions leads to a finite conductivity, a e \ ~ T / ct\n{l / a) , 
at temperatures well below T c . At temperatures above, W ± charge-exchange processes - 
analogous to color exchange through gluons in QCD - effectively stop left-handed charged 
leptons. However, right-handed leptons can carry current, resulting in a e i/T being only a 
factor ~ cos 4 8w smaller than at temperatures below T c . 



I. INTRODUCTION 



The strong magnetic fields measured in many spiral galaxies, B ~ 2 x G Q, are conjectured to 

be produced primordially; proposed mechanisms include fluctuations during an inflationary universe || or 
at the GUT scale Q , and plasma turbulence during the electroweak transition [^|-|| or in the quark-gluon 
hadronization transition . The production and later diffusion of magnetic fields depends crucially on the 
electrical conductivity, a e i, of the matter in the universe; typically, over the age of the universe, t, fields on 
length scales smaller than L ~ (t/iiraei) 1 / 2 are damped. In this paper we calculate a e i in detail below and 
above the electroweak transition scale. 

The electrical conductivity was estimated in [p| in the relaxation time approximation as a e i ~ naT e i/m 
with m ~ T and relaxation time r e i ~ l/(a 2 T), where a = e 2 /4tt. In Refs. and pi the relaxation time 
was corrected with the Coulomb logarithm. A deeper understanding of the screening properties in QED and 
QCD plasmas in recent years has made it possible to calculate a number of transport coefficients including 
viscosities, diffusion coefficients, momentum stopping times, etc., exactly in the weak coupling limit [pd|-p^[; 
also flM. However, calculation of processes that are sensitive to very singular forward scatterings remain 
problematic. For example, the calculated color diffusion and conductivity jig ], even with dynamical screening 
included, remain infrared divergent due to color exchange in forward scatterings. Also the quark and gluon 
damping rates at non-zero momenta calculated by resumming ring diagrams exhibit infrared divergences 
jt6| whose resolution requires more careful analysis including higher order diagrams as, e.g., in the Bloch- 
Nordseick calculation of Ref. jlTj . Charge exchanges through exchange, processes similar to gluon color 
exchange in QCD, are important in forward scatterings at temperatures above the W mass, Mw- While, 
as we show, such processes lead negligible transport by left-handed charged leptons, they do not prevent 
transport by right-handed charged leptons. As a consequence, electrical conduction at temperatures above 
the electroweak transition is large, and does not inhibit generation of magnetic fields JIo[| ; the observed 
magnetic fields in galaxies could then be generated earlier than the electroweak transition j4|J|] and have 
survived until today. More generally, we find that the electrical conductivity is sufficiently large that it does 
not lead to destruction of large-scale magnetic flux over timescales of the expansion of the universe. 



In Sec. II we calculate the electrical conductivity for T <C Myy- In this regime the dominant momentum 
transfer processes in electrical transport arc clcctrodynamic, and one can ignore weak interactions. We turn 
then in Sec. Ill to the very early universe, T ^> T c , where the are effectively massless and their effects 
on electrical conduction must be taken into account. 

II. ELECTRICAL CONDUCTIVITIES IN HIGH TEMPERATURE QED 

We first calculate the electrical conductivity in the electroweak symmetry-broken phase at temperatures 
well below the electroweak boson mass scale, T <C Mw- As we argue below, charged leptons I — e~, t~ 
and anti- leptons i = e + ,/i + ,r + dominate the currents in the regime in which me <C T. In the broken- 
symmetry phase weak interactions between charged particles, which are generally smaller by a factor 
~ (T/Mw) 4 compared with photon-exchange processes, can be ignored. The primary effect of strong in- 
teractions is to limit the drift of strongly interacting particles, and we need consider only electromagnetic 
interactions between charged leptons and quarks. 

Transport processes are most simply described by the Boltzmann kinetic equation for the distribution 
functions, ?ij(p, r, i), of particle species i, of charge e, 

{§t + Vl ' Vr + eE ' Vl ) ni = ~ 2lTV2 |^ l M i2-3 4 | 2 Kn 2 (I ± n 8 )(l ± n 4 ) - n 3 n 4 (l ± m)(l ± n 2 )] 

x (5 Pl +p 2 ,p 3 +p 4 <5(ei + £2 - £3 - £4), (1) 

where E is the external electric field driving the charges, and the right side describes the two-particle collisions 
12 <-> 34 slowing them down. The ± signs refer to bosons and fermions. The sums are over momenta p 2 , 
P3, and P4, and the statistical factor v 2 accounts for the number of leptons and spin projections that 
scatter with particle 1. Massless lepton-lepton or antilepton-antilepton scattering conserves the electrical 
current, and affects the conductivity only in higher order. In lowest order we need to take into account only 
lepton-antilepton scattering. The square of the matrix element for scattering of a charged lepton from initial 
momentum pi to final momentum p 3 with an antilepton from p 2 to p 4 is 

|M 12 ^ 34 | 2 = 2e 4 (s 2 + U 2 )A 2 /(i6eie 2 e 3 e4), (2) 

where s,t,u are the usual Mandelstam variables formed from the incoming and outgoing momenta. 

To solve the kinetic equation for a weak driving field, we write n$ = + $i, where nf — (exp(ej/T) ± 1) _1 
is the global equilibrium distribution (chemical potentials are taken to vanish in the early universe), and 
linearize Eq. ([!]) in the deviations $i of the distribution functions from equilibrium (see Refs. |18[ ], Jl9| and 
Ji"3"| for details); thus Eq. (0) reduces to 

eE- Vl — - = -2TT1S2 J2 |Mi2^34| 2 n?^(l ± ng)(l ± n°) ($1 + $ 2 - $ 3 - $ 4 ) 5 Pl+ p 2 ,p 3+P4 <5(ei + e 2 - e 3 - e 4 ). 
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We set the stage by considering the simple situation of conduction by massless charged leptons (component 
1) and antileptons (component 2) limited by their scattering together; later we include effects of quarks. The 
electric field drives the charged leptons and antileptons into flow with steady but opposite fluid velocities 



Ui = — U2 in the center-of-mass system (see Fig. 1). Assuming that collisions keep the drifting components 
in local thermodynamic equilibrium, we approximate the quasiparticle distribution functions by those of 
relativistic particles in equilibrium with a local flow velocity u, 

"i(Pi) = FT \ 7rp] -r- i i ( 4 ) 

exp[(ej - Ui ■ Pi)/T\ =F 1 

the deviation $i is thus 

$i = -Uj-Piv- ■ 5) 

The ansatz for the distribution function is an excellent approximation; a more elaborate variational 
calculation which we have carried out (see jl3| for analogous calculations for the viscosity) gives almost the 
same distribution function, and a corresponding electrical conductivity only 1.5% smaller than our result, 
Eq. (^) below. Equation (JsJ) gives the total electric current from lepton-antilepton pairs with Nt effectively 
massless species present at temperature T: 

lu = e^chUi = eui^^T 3 , (6) 

where n c h is the number density of electrically charged leptons plus antileptons. Note that since photon 
exchange does not transfer charge (as do W^), particle 3 in Eq. (Q) has the same charge as particle 1, and 
particle 4 the same as 2, and 113 = Ui and 114 = U2 = — ui. 

To calculate a e i we multiply Eq. (^|) by V1V1 and sum over pi, to find an equation of the form: 

- eEN e T 2 /18 = (7) 

where £ results from the right side of ([}]). Since QED interactions are dominated by very singular forward 
scattering arising from massless photon exchange, the sum on the right side in £ diverges unless we include 
Debye screening of longitudinal interactions and dynamical screening of transverse interactions due to Landau 
damping. This is done by including (in the Coulomb gauge) the Coulomb self-energy, I±£, and the transverse 
photon self-energy, Ht, in the longitudinal and transverse electromagnetic propagators. The inverse of the 
longitudinal propagator becomes q 2 + I±£(a;,q), and the inverse of the transverse propagator, t = lu 2 — q 2 
becomes u> 2 — q 2 — Ht(uj, q), where lo and q are the energy and momentum transferred by the electromagnetic 
field in the scattering. The quantity £ can be calculated to leading logarithmic order in the coupling constant 
by expanding for small momentum transfers (see Ref. |l3( ] for details). Small momentum transfer processes 
are screened by 11^ ~ q^ = AiraNiT 2 /3, and TIt ~ iitq 2 D uj / Aq. (Large momentum transfers, q^(p) ~ 3T, 
are cut off by the distribution functions.) The resulting integrals, J q 2 d 2 q/\q 2 + IIz,,7"| 2 , give characteristic 
logarithms, ln(T/g£>), and we find, 

£=^A^ln(C7aiV ; )T 4 . (8) 
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The constant C ~ 1 in the logarithm gives the next to leading order terms (see 113] for the calculation of 
second-order contributions to the viscosity). The electrical conductivity for charged leptons is thus p2] , [l9[ | 

•S'-^-^SE^j. ».«r«r aCD . (9) 



Note that the number of lepton species drops out except in the logarithm. The above calculation taking 
only electrons as massless leptons (Ni = 1) gives a first approximation to the electrical conductivity in the 
temperature range m e <C T <C Tqgp-, below the hadronization transition, Tqgp ~ 150 GeV, at which 
hadronic matter undergoes a transition to a quark-gluon plasma. Thermal pions and muons in fact also 
reduce the conductivity by scattering electrons, but they do not become significant current carriers because 
their masses are close to Tqgp- 

For temperatures T > Tqgp, the matter consists of leptons and deconfined quarks. The quarks themselves 
contribute very little to the current, since strong interactions limit their drift velocity. The quark drift velocity 
can be estimated by replacing a by the strong interaction fine structure constant, a s , in £, Eqs. (]?]) and (^), 
which yields u q ~ u^(a 2 In a -1 ) /{a 2 s lna^ 1 ). 

Even though quarks do not contribute significantly to the currents they are effective scatterers, and thus 
modify the conductivity (an effect ignored in the recent numerical analysis of Ref. j2(J). To calculate the 
quark contribution to the lepton conductivity, we note that the collision term between leptons (1,3) and 
quarks (2,4) includes the following additional factors compared with lepton-lepton scattering: a factor 1/2, 
because the quark velocity, 112, is essentially zero, a factor 3 from colors, and a factor 2 because the charged 
leptons collide on both q and q; finally we must sum over flavors with a weight Q 2 q , where Q q e is the charge 
of quark flavor q — u,d, s, c, b, t. We must also divide by the number of leptons, JVj, to take into account the 
number of quark scatterings relative to lepton scatterings. Including Iq and tq collisions on the right side of 
Eq. @ we find the total electrical conductivity of the early universe Jl2|,[l9| : 

n Nl Nl 3C(3) T r *t*m rim 

el ~n 1 + 3 E? Q\ * ~n 1 + 3 Q\ ™T « Hihm ' t ^ P «t«m w . (10) 

The charged lepton and quark numbers N[ and N q count only the species present in the plasma at a given 
temperature, i.e., those with masses rrii^T. Figure 2 illustrates the conductivities (p||lC|). For simplicity this 
figure assumes that the quarks and leptons make their appearance abruptly when T becomes > mf, in reality 
they are gradually produced thermally as T approaches their masses. Since a range of particle energies in 
the neighborhood of the temperature is included, possible resonance effects in scatterings are smeared out. 

We will not attempt to calculate the electrical conductivity in the range Mw^T^JT c below the critical 



temperature. Recent lattice calculations 21 predict a relatively sharp transition at a temperature T c ~ 100 
GeV from the symmetry broken phase at low temperatures, T -C T c , to the symmetric phase at T > T c . 
The transition is sensitive, however, to the unknown Higgs mass, with a first order transition predicted only 
for Higgs masses below ~ 90GeV. The calculations of the conductivity are technically more involved when 
masses are comparable to the temperature. Furthermore one must begin to include contributions of the 
to currents and scatterings, as the thermal suppression of their density decreases near the transition. 

III. THE SYMMETRY-RESTORED PHASE 



To calculate the conductivity well above the electroweak transition, T ^> T c , where the electroweak 
symmetries are fully restored, we describe the electroweak interactions by the standard model Weinberg- 
Salam Lagrangian with minimal Higgs: 
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H 2 (f^(f> - AO^) 2 - G x L^ii + iG 2 LT 2 (j)*R + h.c. . (11) 



Here L denotes left-handed doublet and R right-handed singlet leptons or quarks, e = gsmOw — g' cosQw 
and electrical charge Q = T3 + y/2. The last terms provide masses for leptons and quarks in the low 
temperature phase, T < T c , where the Higgs field has a non-vanishing vacuum expectation value {4>) = 
(0, v)/V2; at zero temperature v 2 = -fi 2 /X = 1/(G F V2) = AM^/g 2 = (246GeV) 2 . 

At temperatures below T c , where fi 2 < 0, the Higgs mechanism naturally selects the representation W^, 
Z° , and 7 of the four intermediate vector bosons. At temperatures above the transition - where (cj)) vanishes 
for a sharp transition, or tends to zero for a crossover - we consider driving the system with external vector 
potentials A a — B, W , W 3 , which give rise to corresponding "electric" fields E Q , where 

E? = Ft = diA% - d Af, A a = B (12) 
= F? = diA% - d A° - ge abc A\A c , A a = W l ,W 2 ,W\ (13) 

One can equivalently drive the system with the electromagnetic and weak fields derived from A, Z°, and 
W ± , as when T < T c , or any other rotated combination of these. We consider here only the weak field limit 
and ignore the nonlinear driving terms in Eq. (fi"3"|). The self-couplings between gauge bosons are important, 
however, in the scattering processes in the plasma determining the conductivity, as we discuss below j^] 

The electroweak fields A b act on the matter to generate currents J a of the various particles present in the 
plasma, such as left and right-handed leptons and their antiparticles, and quarks, vector bosons, and Higgs 
bosons. The Higgs and vector boson contributions are, as we shall see, negligible. Therefore the significant 
terms in the currents are 

j b = i^iL^YL + Ry^YR) (14) 

J&i = \^nh. (15) 

We define the conductivity tensor a a b in general by 

J a = a ab E b . (16) 

Equation (|l^) with the equations of motion in a gauge with d^A 1 ^ = yields the weak field flux diffusion 
equation for the transverse component of the fields, as in QED, 



1 Pure SU(2) gauge fields undergo a confining first order phase transition [g2| at a critical temperature, T c , where the 
fields are intrinsically strongly interacting. In an academic universe without Higgs fields the electroweak phase transi- 
tion is non-existent. If, in this case, we run the temperature from far above Mw to low temperatures the electroweak 
running coupling constants diverge at a very small temperature ~ 10 -23 GeV, signalling that the interactions of the 
fields have become strong, and that one is in the neighborhood of the SU(2) phase transition. For our purposes we 
are therefore safe in ignoring such non-linear effects and the SU(2) phase transition. However, the nature of electrical 
conduction in the confined state by charged "mesons" formed, e.g., by e + v e or ud, remains an interesting problem in 
principle. We thank Eduardo Fradkin for calling this issue to our attention. 



(d 2 - V 2 )A a = a ab d t A b . 



(17) 



describing the decay of weak fields in terms of the the conductivity. 

The electroweak U(l) x SU (2) symmetry implies that the conductivity tensor, a ab , in the high temperature 
phase is diagonal in the representation a, b = B, W 1 , W 2 , W 3 , as can be seen directly from the (weak held) 
Kubo formula 

(Jab = - lhn lim — Im ( J a J & ) irr , (18) 

which relates the conductivity to (one-boson irreducible) current-current correlation functions.^ The con- 
struction of the conductivity in terms of the Kubo formula assures that the conductivity and hence the 
related entropy production in electrical conduction are positive. Then 

/ (tbb 
a ww 
a ~ aww 
V aww 

Due to isospin symmetry of the VF-interactions the conductivities aw i w i are the same, = aww, but differ 
from the B-held conductivity, asB- 

The calculation of the conductivities asB and aww m the weak held limit parallels that done for T « T c . 
The main difference is that weak interactions are no longer suppressed by a factor (T/Mw) 4 and the exchange 
of electroweak vector bosons must be included. The conductivity, agg, for the abelian gauge held B can 
be calculated similarly to the electrical conductivity at T <C T c . Taking into account the fact that both 
left-handed neutrinos and charged leptons couple to the -B-field with the same sign, and that they scatter 
the same way, their flow velocities are equal. Consequently, in the scatterings 12 <-> 34, ui = 113 and U2 = U4, 
whether or not the interaction is by charge exchange. The situation is thus similar to electrodynamic case. 

Although the quarks and are charged, their drifts in the presence of an electric held do not signihcantly 
contribute to the electrical conductivity. Charge flow of the quarks is stopped by strong interactions, while 
similarly flows of the are effectively stopped by W + + W~ — > Z°, via the triple boson coupling. Charged 
Higgs bosons are likewise stopped via W^^cj) couplings. These particles do, however, affect the conductivity 
by scattering leptons. 

The lepton and quark mass terms in the Weinberg-Salam Lagrangian provide masses only when the Higgs 
held has a non-zero expectation value. For T » T c the quarks and leptons have thermal masses, which, 
for the longitudinal (electric) degrees of freedom, are of order the plasma frequency, m p i ~ gT, and of 
likely order m mag ~ g 2 T |23f| for the transverse (magnetic) mass. These small masses give rise to spin-flip 
interactions changing the helicity; such interactions are, however, suppressed by factors of m/T, and can 
therefore be neglected here. The mass terms also provide a small coupling, Gi = \f2raijv, between the Higgs 



2 Since the linearized currents are proportional to the Pauli spin matrices r a for the W (a=l,2,3) fields and the 
identity matrix to = 1 for the B field, one finds explicitly in one-loop order that a a b oc Tr{r a Ti,} = 28 a b (a=0,l,2,3). 
Including the dominant interactions in the current-current loop is basically equivalent to solving the Boltzmann 
equation, which produces no off-diagonal elements. 




and leptons, proportional to the ratio of the lepton mass to the vacuum expectation value, which leads to a 
negligibly small contribution to the conductivity. Even the coupling of r mesons with the Higgs is a factor 
G 2 /e 2 ~ 10 -3 smaller than their coupling to the B field. Charge transfer scatterings via Higgs exchange 
are more singular than scatterings via B exchange, and are enhanced by a factor 1/e 2 ; nonetheless such 
processes are negligible. 

These considerations imply that the B current consists primarily of right-handed , fi^ and , interact- 
ing only through exchange of uncharged vector bosons B, or equivalently 7 and Z°. Because the left-handed 
leptons interact through W as well as through B, they give only a minor contribution to the current. They 
are, however, effective scattcrers of right-handed leptons. The resulting conductivity is 

obb = \ jWfo a f ) = 9 cos2 9w a^ , T » T c , (20) 

where the Y qa L arc the right and left-handed quark hypercharges, and the Yr,/, are the right and left-handed 
charged lepton hypercharges. The terms entering the prefactor are: i) a factor 1/2 because only the right- 
handed leptons contribute significantly to the conductivity; ii) a net factor of cos 2 9w because the B field 
coupling to right-handed leptons and the current Jb contain factors of g' = e/ cos 9wj while the square of the 
matrix element contains a factor (e/ cos9w) 4 ', and iii) a factor (Yj + 2Y£)/8 = 3/4 in the scatterings of the 
right-handed charged leptons with right and left-handed leptons, and iv) a factor ^2^ q (Y 2 R ~^~Yq L )/& = H/12 
in the scatterings of the right-handed charged leptons with right and left-handed quarks. The factor 9/19 
holds in the limit that the temperature is much greater than the top quark mass, mf, excluding the top 
quark for T c < T < rrit gives 108/211 instead. 

Applying a W 3 field to the electroweak plasma drives the charged leptons and neutrinos oppositely since 
they couple through 5T3W3. In this case, exchanges of W dominate the interactions as charge is transferred 
in the singular forward scatterings, so that 113 = 112 = — Ui. The collision term is then weighted by a factor 
(P1+P2) instead of a factor (pi— P2) = q and one ends up with an integral J p 2 dq 2 / (q 2 +q' D ,) 2 ~ T 2 jq 2 u ~ a -1 
for the longitudinal part of the interaction. For the transverse part of the interaction one encounters a 
logarithmic singularity; while Landau damping is not sufficient to screen the interaction, a magnetic mass, 
fnmag-, wm provide an infrared cutoff. Besides the logarithms, the factor oT 1 remains and we expect that 

cr\vw ~ uobb- (21) 

This effect of H /± exchange is analogous to the way gluon exchange in QCD gives strong stopping and 
reduces the "color conductivity" significantly p5[ ; similar effects are seen in spin diffusion in Fermi liquids 
& 

The electrical conductivity is found from obb and <Jww by rotating the B and W 3 fields and currents by 



the Weinberg angle; using Eq. (16) we obtain 



) = K(6 w )aK(-6 w ) ( ^ 



a B B cos 2 8 W + a ww sin 2 9 W {&bb &ww) cos 6\y sin 6\y \ ( A 
[pww — &bb) cos 9w sin 8\y ct B b sin 2 9w + oww cos 2 9w J \ % 



(22) 



Thus the electrical conductivity is given by 



(J A A — obb cos 2 6 w + &ww sin 2 B w \ 



(23) 



<J e i/T above the electroweak transition differs from that below mainly by a factor ~ cos 4 9w ~ 0.6. 

In the wide temperature range we are considering the coupling constants in fact run as <x;(Q) = on([x) + 
bi\n(Q//j,) where the coefficients hi are found by renormalization group calculations [^4| , ^5"1| . In a high 
temperature plasma typical momentum transfers Q are of order qu ~ eT. The exact values employed for Q 
is not important as the couplings only increase logarithmically with temperature. In the temperature range 
1 to 10 6 GeV, a" 1 varies from 130 to 123 and sin 2 6 W from 0.21 to 0.28. 



We have calculated the electrical and electroweak conductivities in the early universe over a wide range 
of temperatures. Typically, a e i ~ T/a 2 ln(l/cc), where the logarithmic dependence on the coupling constant 
arises from Debye and dynamical screening of small momentum-transfer interactions. In the quark-gluon 
plasma, at T ^> Tqgp ~ 150 MeV, the additional stopping on quarks reduces the electrical conductivity 
from that in the hadronic phase. In the electroweak symmetry-restored phase, r>T c , interactions between 
leptons and W ± and Z° bosons reduce the conductivity further. The electrical conductivity does not vanish 
(as one might have imagined to result from singular unscreened VT^-exchanges), and is larger than previous 
estimates, within an order of magnitude. The current is carried mainly by right-handed leptons since they 
interact only through exchange of 7 and Z . 

From the above analysis we can infer the qualitative behavior of other transport coefficients. The charac- 
teristic electrical relaxation time, r e / ~ (or lxi(l / a)T) , defined from a ~ e 2 nr e ;/T, is a typical "transport 
time" which determines relaxation of transport processes when charges are involved. Right-handed leptons 
interact through Z a exchanges only, whereas left-handed leptons may change into neutrinos by W exchanges 
as well. Since Z° exchange is similar to photon exchange when r>T c , the characteristic relaxation time is 
similar to that for electrical conduction, t v ~ (a 2 ln(l/a)T) _1 (except for the dependence on the Weinberg 
angle). Thus the viscosity is 77 ~ t v ~ T 3 /(a 2 ln(l/a)). For T <gC Mw the neutrino interaction is suppressed 
by a factor (T /My/) 41 ', in this regime neutrinos have longest mean free paths and dominate the viscosity. Jl£| ] 

The electrical conductivity of the plasma in the early universe is sufficiently large that large-scale magnetic 
flux present in this period does not diffuse significantly over timescales of the expansion of the universe. The 
time for magnetic flux to diffuse on a distance scale L is t<#// ~ a e iL 2 . Since the expansion timescale t exp 
is ~ l/(ipianck^ 2 ), where tpianck ~ 10 -43 s is the Planck time, one readily finds that 



where x = L/ct exp is the diffusion length scale in units of the distance to the horizon. As described in Refs. 
fiof and ||, sufficiently large domains with magnetic fields in the early universe would survive to produce 
the primordial magnetic fields observed today. 
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IV. SUMMARY AND OUTLOOK 
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Fig. 1: The particle currents generated by an electric field in the early universe at T ^> T c . Right-handed 
charged leptons, I — e, /it, interact only through 7 and Z° exchanges while left-handed leptons also interact 
by W ± exchanges, which drag neutrinos along and decrease the current. The vector and Higgs bosons, 
7, Z°, , <fi, cannot flow due to W ± exchanges, and qqg are stopped by strong interactions. 

Fig. 2: Electrical conductivity vs. temperature. The temperatures, where the transitions from hadronic to 
quark-gluon plasma and electroweak symmetry breaking occur, are indicated by QGP and EW respectively. 
The conductivity a e i is given by Eqs. (|9|,[l0|,^0|) in the three regions and are extrapolated into the regions of 
the phase transitions. The quark and lepton masses in the figure indicate the temperatures at which they 
are thermally produced and thus affect the conductivity (see Eq. ( |Io| ) and discussion in text). 
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